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Jj 要 : 过 去 70a 我 国 在 干旱 半 干 旱地 区 展开 大 规模 植树 造林 活动 并 取得 显著 成 效 ,由 于 受 人 类 活动 和 气候 影 
响 ,导致 局 部 地 区 造林 失败 。 虽 然 在 干旱 地 区 造林 的 洪 力 被 得 到 认可 ,但 是 造林 具体 应 集中 在 哪些 微观 地 区 尚未 
可 知 。 本 文 以 黄土 高 原 典型 区 域 山西 省 为 例 ,构建 宜 林 地 识别 框架 ,利用 PLUS 模型 和 马尔 可 夫 链 模拟 2030 年 不 同 
情景 下 土地 利用 变化 ,分析 林 地 扩张 、 林 分 内 部 变化 趋势 和 空间 分 布 。 结 果 表 明 : 当前 山西 省 林地 增长 空间 为 
5.38% ,上 且 东南 部 为 主要 潜在 增长 区 ,中 西部 林地 边缘 存在 有 林地 退化 的 可 能 ,而 北部 林地 破碎 化 程度 较 高 ; 林 分 内 
部 互相 转换 频率 较 高 ; 按 扩张 潜力 ,依次 为 有 林地 > 灌木 林 > 路 林地 > 其 他 林地 。 本 研究 可 为 山西 省 植树 造林 有 效 管 


理 和 森林 固 碳水 平 提升 提供 参考 。 
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过 去 60 a 全 球 干旱 半 干 旱地 区 土地 退化 面积 不 
汤 扩大 ,平均 每 年 退化 面积 约 有 1000 hm?" eh BE 
的 人 类 活动 以 及 自然 因素 导致 干旱 半 干 旱地 区 土 
地 退化 严重 ,当前 退化 土地 的 恢复 成 为 干旱 半 干 旱 
地 区 面临 的 重大 挑战 之 一 。 在 此 背景 下 ,许多 国家 
展开 了 包括 植树 造林 在 内 的 森林 恢复 项 目 , 旨 在 组 
解 因 干 旱 而 造成 的 地 方 经 济 损失 ”> 。 

森林 恢复 是 维护 碳 平衡、 缓解 全 球 气候 变 暖 的 
有 效 策 略 之 一 。 森 林 生 态 系统 在 加 强 碳 储量 .提升 
生态 系统 承载 力 ,维护 生态 系统 平衡 等 方面 发 挥 着 
重要 作用 *"。 由 于 土地 过 度 利用 、 林 业 无 序 开发 、 
自然 灾害 等 原因 ,导致 全 球 森 林 面 积 净 损 失 达 80x 
10 km 。 有 研究 表明 ,在 当前 气候 条 件 背 景 下 ,全 


综 上 ,围绕 干旱 半 干 旱地 区 造林 的 重要 课题 ， 
已 有 研究 主要 从 以 下 几 个 方面 进行 了 论述 :第 一 个 
核心 议题 是 关于 造林 的 适宜 性 。 对 于 在 植树 造林 
过 程 中 的 适宜 度 ,学 术 界 主要 有 以 下 两 种 观点 。 一 
是 植被 选择 的 适宜 性 … ,在 干旱 半 干 旱地 区 ,外 来 
单一 植被 的 种 植 往往 会 导致 后 期 为 维护 生物 多 样 
性 而 产生 高 兄 成 本 ,造成 不 必要 的 政府 负担 ”。 二 
是 造林 规模 的 适宜 性 ,大 规模 盲目 的 森林 扩张 可 能 
更 易 受 到 本 地 干旱 化 的 影响 ,包括 在 不 适宜 的 地 点 
进行 大 规模 的 造林 ”。 以 上 两 种 观点 表明 不 适宜 
的 人 工 修复 会 造成 后 期 森林 管 护 成 本 增加 。 虽 然 
学 术 界 越 来 越 意识 到 这 一 问题 的 严重 性 ,但 是 目前 
国内 关于 在 微观 地 区 造林 的 空间 数据 和 信息 尚 不 


球 约 有 9x10’ km 的 森林 增长 空间 ”。 其 中 ,在 干旱 
半 干 旱地 区 森林 增长 的 潜力 巨大 , 约 占 潜在 森林 恢 
复 面 积 的 36%~42%™ ,但 这 可 能 是 相对 乐观 的 估 
计 。 实 际 上 ,在 干旱 半 干 旱地 区 造林 面临 种 种 困 
难 , 虽 然 许 多 国家 和 组 织 积极 倡导 森林 恢复 的 重要 
性 ,但 是 面 对 当 前 紧迫 的 全 球 变 暖 趋势 ,适宜 干旱 
半 干 旱地 区 森林 恢复 的 潜在 位 置 在 哪 ,以 及 如 何 恢 
复 等 问题 还 没有 确切 的 定论 ""。 
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完善 ,这 导致 无 法 准确 掌握 干旱 半 干 旱地 区 造林 的 
空间 适宜 与 否 。 

第 二 个 核心 议题 是 关于 干旱 半 干 旱地 区 造林 
的 贡献 度 。 在 诸多 学 术 人 研究 中 ,关于 旱地 森林 的 人 研 
究 往 往 少 于 其 他 生物 群落 ,虽然 与 其 他 生物 群落 相 
比 旱 地 和 森林 的 植被 覆盖 率 极 低 ,但 是 其 在 全 球 碳 循 
环 中 发 挥 着 重要 作用 ,部 分 旱地 和 森林 的 碳 密 度 可 能 
高 于 或 近 于 临近 森林 生物 群落 的 碳 密度 。 我 国 
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黄土 高 原 即 为 典型 代表 区 域 , 该 地 区 森林 植被 的 恢 
复 显著 提升 了 区 域 总 碳 储量 ”。 为 此 ,我 国 根据 
《中 国 森 林 经 营 规划 (2016 一 2050)》(http://www.for- 
estry.gov.cn/) , 需 进一步 在 干旱 半 干 旱地 区 建立 大 量 
人 工 林 来 实现 碳 中 和 目标 。 而 在 干旱 区 域 实施 造 
林 工 程 ,需要 考虑 到 因 造 林 而 涉及 在 生态 ,经济 和 
社会 复合 生态 系统 中 的 权衡 问题 。 因 此 ,潜在 造林 
面积 的 调查 和 造林 区 域 的 识别 对 干旱 半 干 旱地 区 
造林 具有 重大 意义 。 

第 三 方面 需要 重点 关注 我 国 对 植树 造林 的 政 
策 把 控 。 和 森林 恢复 被 视 为 是 国家 自主 贡献 的 主要 
组 成 部 分 。 从 我 国 第 三 次 全 国土 地 调查 来 看 ,我 国 
政府 对 当前 的 生态 修复 似乎 越 来 越 保持 谨慎 的 态 
度 ,严格 把 控 耕 地 向 林地 .草地 和 园地 的 转换 成 为 
主要 趋势 ,退耕 还 林 的 表述 也 在 进一步 被 弱化 , 反 
而 更 加 强调 耕地 保护 红线 的 重要 性 。 因 此 ,对 于 干 
旱 半 干旱 地 区 造林 的 实现 更 加 要 求 对 区 域 选 择 的 
重要 性 和 造林 效率 的 提升 。 

综 上 ,梳理 以 上 三 方面 的 研究 对 我 国 干旱 半 干 
旱地 区 造林 具有 重要 意义 。 黄 土 高 原作 为 我 国 干 
旱 半 干旱 地 区 典型 的 造林 区 域 ,其 森林 生态 修复 受 
到 社会 各 界 关 注 。 目 前 ,关于 黄土 高 原 地 区 植树 造 
林 的 研究 主要 集中 在 关于 造林 过 程 中 的 变化 1、 
造林 的 相关 影响 因素 评析 ”以 及 造林 成 效 E 
方面 。 少 部 分 研究 集中 在 以 气候 条 件 和 政策 为 背 
景 ,对 造林 潜力 评估 人 研究。 但 仍 面临 着 具体 方法 
论 方面 的 困扰 :首先 ,针对 森林 生态 修复 的 研究 很 
少 体 现在 森林 恢复 过 程 中 关于 林 分 内 部 变化 , 例 
如 , 某 些 区 域 受 气候 和 人 类 活动 影响 ,森林 内 部 出 
现 灌木 林 向 琉 林地 退化 的 趋势 ,而 这 在 森林 总 体面 
积 上 却 并 未 表现 出 来 ,但 显然 森林 质量 是 退化 的 ; 
其 次 ,相关 研究 关于 受 人 类 活动 干扰 导致 土地 利用 
变化 的 森林 动态 研究 也 较 少 体现 。 因 此 ,结合 上 述 
几 个 方面 关于 造林 的 相关 论证 ,本 文 认为 , 因 人 为 
因素 促使 土地 利用 变化 下 的 森林 生态 恢复 的 具体 
造林 潜力 区 还 需 进一步 挖掘 ,在 森林 恢复 过 程 中 关 
于 林 分 内 部 的 变化 也 需 进 一 步 分 析 。 

山西 省 是 黄土 高 原 地 区 人 工 造林 典型 区 域 , 境 
内 宜 林地 面积 广阔 。 近 年 来 ,在 耕地 保护 政策 下 耕 
地 向 林地 .草地 转化 面积 的 减少 限制 了 退耕 还 林 
〈 草 ) 的 发 展 空间 。 此 外 ,自然 因素 和 人 为 管理 不 善 
也 导致 林 分 质量 下 降 。 因 此 ,本 研究 以 山西 省 为 主 


要 研究 区 域 , 运 用 PLUS 模型 和 马尔 科 夫 链 对 不 同 
情景 下 林地 利用 结构 变化 进行 模拟 ,特别 对 林 分 内 
部 之 间 转 化 进行 细 分 。 本 研究 目的 :(1) 评 佑 林地 
生长 的 区 域 适 宜 程 度 ;(2) 评估 山西 省 林地 增长 空 
间 及 其 分 布 情况 ;(3) 识别 林地 优先 发 展区 域 。 


1 研究 区 域 和 数据 来 源 


1.1 研究 区 概况 

山西 省 位 于 我 国 中 部 地 区 (34?34' ~40°44'N, 
110°14’~114°33’E) ,地 处 第 二 阶梯 , 东 起 太行 山 , 西 
至 吕梁 山 ,中 部 为 汾 河 河 谷 ,形成 “两 山 夹 一 川 ” 的 
地 形 特点 。 地 理 单 元 相对 闭塞 ,地 势 由 东北 逐渐 向 
西南 方向 下 降 , 主要 景观 单元 为 山地 高 原 、 山 地、 丘 
陵 。 为 温带 大 陆 性 季风 气候 ,全 年 降水 季节 分 配 不 
均 , 介 于 358~621 mm ,水 资源 相对 较 少 。 省 域内 的 
主要 生态 系统 为 森林 生态 系统 .湿地 生态 系统 以 及 
草原 生态 系统 ,但 生态 系统 总 体 较为 脆弱 ” 。 
1.2 数据 来 源 

本 研究 数据 主要 包括 土地 利用 数据 和 土地 利 
用 变化 驱动 因子 。 其 中 ,土地 利用 类 型 数据 来 自 中 
国 科学 院 资 源 环 境 科 学 与 数据 中 心 ( 分 辩 率 30 mx 
30 m)。 为 直观 分 析 山 西 省 林地 内 部 转换 情况 ,各 
土地 利用 类 型 为 有 林地 灌木 地 Lit PRO. 、 其 他 林 
地 耕地 、 草 地、 水域. 建设 用 地 以 及 未 利用 地 。 对 
于 宜 林 地 识别 的 土地 利用 需求 驱动 因子 主要 分 为 3 
大 类 ( 表 1), 即 土壤 驱动 因子 (9 项 ) 、 社 会 经 济 驱 动 
因子 (10 项 ) 以 及 气候 环境 驱动 因子 (5 项 )。 其 中 ， 
坡度 和 坡 向 数据 根据 DEM 数据 计算 获得 ,气温 和 降 
水 数据 主要 采用 1970 一 2000 年 的 均值 代替 。 


2 模型 与 方法 


本 研究 基于 马尔 科 夫 链 .PLUS 模 型 构建 了 一 个 
综合 识别 框架 (图 1), 首 先 ,采用 马尔 科 夫 链 预 测 山 
西 省 不 同情 景 下 土地 需求 预测 值 。 其 次 ,利用 PLUS 
模型 将 预测 值 和 各 土地 利用 驱动 因子 分 配 到 不 同 
土地 利用 类 型 中 ,获得 多 情景 下 土地 利用 空间 分 布 
情况 并 进行 对 比分 析 。 

2.1 基于 土地 利用 需求 数量 和 空间 约束 的 多 情景 
设置 

本 研究 利用 马尔 可 夫 链 由 在 依据 可 持续 发 展 
目标 (SDGs)"™ ,设置 不 同情 景 , 并 对 不 同情 景 下 
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表 1 土地 利用 需求 变化 驱动 因子 


Tab.1 Drivers of land use demand change 


数据 类 型 数据 名 称 年 份 分 辩 率 /m 数据 来 源 

土壤 驱动 因子 土壤 融 水 量 ; 土 壤 PH 值 ; 距 基 黏 岩 深 度 ; 有 机 土壤 的 累 2017 250 hittps://www.isric.org/ 
积 概率 ;土壤 有 机 碳 储量 ; 含 沙 量 ; 含量 ;土壤 纹理 
土壤 类 型 1995 1000 http://www.resde.cn/ 

社会 经 济 驱 动因 子 人 口 ;GDP; 临 近 城 市 ;临近 农村 ; ea a 2015 30 https://www.openstreetmap.org/ 
公路 ;临近 主干 道 ;临近 二 级 公路 ;临近 三 级 公路 ;临近 
9 级 公路 

气候 及 环境 驱动 因子 DEM ;坡度 ; 坡 向 2016 90 https://Ipdaac.usgs.gov/ 
年 平均 温度 ;年 降水 量 1970 一 2000 30 http://www.worldclim.org/ 
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v 
基准 情景 || 生态 保护 情景 | | ”城镇 化 情景 
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2030 年 山西 省 各 情景 
中 不 同 地 类 需求 面积 
2010 年 和 2020 年 实 


PLUS 模型 
际 土地 利用 栅 格 图 
模拟 不 同情 景 下 土 土地 利用 | THAS k— 
地 利用 空间 分 布 扩张 部 分 || 于 地 利用 变化 驱动 | 
i 因子 | 


图 1 PLUS 模型 和 马尔 可 夫 链 模型 集成 框架 
Fig. 1 Integration framework of PLUS model and Markov chain mode 


不 同 土地 利用 类 型 图 
(20104E., 20204F ) 


未 来 土地 利用 需求 面积 进行 预测 和 设置 。 不 同情 。 用 类 型 向 男 一 种 土地 利用 类 型 转变 的 权限 ,“1” 表 
景 下 土地 利用 需求 值 的 变化 主要 通过 改变 2010 一 。 示 可 以 转换 ,“0” 表 示 禁 止 转换 。 以 2030 年 为 目标 
2020 年 马尔 可 夫 链 中 土地 利用 类 型 间 的 转换 概率 。 年 ,设置 3 种 发 展 情景 来 预 佑 未 来 山西 省 不 同 发 展 
实现 ( 表 2)。 此 外 ,转换 约束 矩阵 表明 一 种 土地 利 ”情景 下 土地 利用 变化 。 


表 2 2010 一 2020 年 土地 利用 类 型 转移 概率 矩阵 
Tab.2 Probability matrix of land use type transfer between 2010 and 2020 


1 2 3 4 5 6 7 8 9 
1 0.8888 0.0271 0.0103 0.0016 0.0298 0.0375 0.0010 0.0038 0.0001 
2 0.0277 0.8426 0.0127 0.0009 0.0412 0.0670 0.0017 0.0060 0.0001 
3 0.0327 0.0283 0.7247 0.0018 0.0931 0.1073 0.0012 0.0107 0.0002 
4 0.1043 0.0468 0.0278 0.6230 0.1050 0.0439 0.0063 0.0424 0.0005 
5 0.0134 0.0139 0.0114 0.0051 0.7907 0.1082 0.0055 0.0517 0.0002 
6 0.0234 0.0294 0.0142 0.0016 0.1400 0.7723 0.0020 0.0167 0.0003 
7 0.0165 0.0108 0.0086 0.0056 0.2388 0.0825 0.5666 0.0641 0.0064 
8 0.0028 0.0032 0.0026 0.0017 0.1748 0.0213 0.0039 0.7894 0.0003 
9 0.0084 0.0084 0.02633 0.0098 0.2453 0.1266 0.0620 0.0867 0.4265 


注 : 表 中 1 表示 有 林地 ;2 表示 灌木 林 ;3 GS TRH 5 4 Ro AAR 5 5 CIS EHI 5 6 CS HH 5 7 RIKER; 8 o LED MB 9 Ro AR A FH 


下 同 。 
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(1) 基准 情景 :该 情境 下 ,对 土地 利用 不 设置 约 
R A& d ,假设 2020 一 2030 年 的 土地 转移 概率 和 
2010 一 2020 年 的 转移 概率 相同 ,人 口 及 经 济 发 展 模 
式 继续 延续 当前 状态 发 展 。 该 情景 模拟 了 无 外 部 
约束 力 下 的 土地 利用 发 展 情况 。 

(2) 生态 保护 情景 :按照 《全 国 森 林 经 济 规划 
(2016 一 2050 年 )》 和 《山西 省 国土 空间 规划 (2020 一 
2035 年 )) 的 要 求 设置 生态 保护 情景 ,山西 省 林地 将 
被 进一步 保护 ,对 林地 向 其 他 地 类 的 转换 进行 限制 
( 表 3)。 在 该 情景 中 2020 一 2030 年 4 类 林地 转 入 概 
率 将 比 2010 一 2020 年 高 , 转 出 概率 下 降 , 同 时 限制 
建设 用 地 的 扩张 ( 表 4)。 


(3) 城镇 化 发 展 情景 :在 城镇 化 情景 下 ,建设 用 
地 会 过 度 占 用 耕地 .草地 等 成 为 必然 。 依 据 这 一 实 
际 , 相 比 于 2010 一 2020 年 土地 转移 概率 ,本 研究 将 
各 类 林地 耕地 草地、 未 利用 地 转 和 人 建设 用 地 的 概 
率 提高 至 20% ,建设 用 地 向 其 他 土地 利用 类 型 转 出 
的 概率 降低 30% (排除 耕地 )( 表 5)。 土 地 转换 矩阵 
与 基准 情景 的 土地 转换 和 矩阵 一 致 ( 表 6)。 
22 PLUS 模型 

PLUS 模型 主要 由 两 个 基本 模块 组 成 ,分 别 为 用 
HOP SKA) AT EIS LEAS fI CARS BERL), LEAS $E 
略 通过 提取 两 期 土地 利用 类 型 扩张 部 分 ,借助 随机 
森林 法 对 扩张 部 分 采样 并 分 析 其 驱动 力 ;CARS 模 


表 3 生态 保护 情景 下 各 土地 利用 类 型 转换 矩阵 


Tab.3 Conversion matrix for each land use type under the ecological conservation scenario 


序号 1 2 3 4 5 6 7 8 9 
1 1 1 1 0 0 0 0 0 0 
2 1 1 1 0 0 0 1 0 0 
3 1 1 1 1 0 1 0 0 0 
4 1 1 1 1 1 1 1 1 0 
5 1 1 1 1 1 1 1 1 0 
6 1 1 1 0 0 1 1 1 0 
7 1 1 0 0 0 0 1 0 0 
8 0 0 0 0 0 0 0 1 0 
9 1 1 1 1 1 1 1 1 1 
#4 生态 保护 情景 下 各 土地 利用 类 型 转换 比重 
Tab.4 Share of land use conversion by type under the ecological conservation scenario 
序号 1 2 3 4 5 6 7 8 9 
1 = = = 20% -20% -20% = -30% -30% 
2 20% - - 一 -20% -10% € -30% -30% 
3 20% = = = -10% -10% = -20% -20% 
4 20% e = = -10% -10% = -20% -20% 
5 20% = B - - B 5 - - 
6 m 20% 20% = -20% -20% = = = 
7 = = = = = zm = E = 
8 = = = = = E = = = 
9 10% 20% 20% 20% = 10% 10% = = 
5 城镇 化 发 展 情景 下 各 土地 利用 类 型 转换 比重 
Tab.5 Share of land use conversion by type under the urbanisation development scenario 
号 1 2 3 4 5 6 7 8 9 
= -20% -20% -20% -30% -20% - -10% -30% 
20% - 一 -20% -20% -20% - 一 -20% 
20% = -20% -20% -20% = = -20% 
= = 20% i = i = 20% 一 
1096 20% 20% 20% = = = 20% = 


- 20% 20% - 
-30% -30% 
20% 20% 
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表 6 城镇 化 发 展 情景 下 各 土地 利用 类 型 转换 矩阵 


Tab. 6 Conversion matrix for each land use type under urbanization development scenarios 
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SEF BENLER EAE SCA ER OL i), AC EI E 
系数 .领域 效应 和 开发 概率 来 模拟 未 来 土地 利用 

需要 注意 的 是 ,在 本 研究 中 PLUS 模型 的 相关 
参数 设置 ,LEAS 模 块 中 ,随机 采样 为 50, 采 样 率 为 
0.01。 在 CARS 模 型 中 , 斑 块 生成 阔 值 为 0.9 ,随机 斑 
块 种 子 概率 为 0.1, 随 机 种 子 比例 为 0.001, 邻 域 范 围 
为 3。 


在 模拟 之 前 , 须 利 用 PLUS 模 型 对 山西 省 土地 
利用 模拟 进行 精度 验证 。 首 先 ,选取 有 影响 土地 利用 / 
土地 覆盖 变化 的 26 种 因素 输入 到 模型 中 ,并 生成 9 
种 土地 利用 类 型 的 适宜 性 概率 (图 2)。 

2.3 模型 精度 验证 

为 保证 模拟 的 精准 性 ,利用 Kappa 系数 对 PLUS 
模型 的 可 靠 性 进行 确定 , A kappa 系数 介 于 0.6~0.8， 
则 表明 仿真 结果 较 高 , 若 介 于 0.8~1.0, 则 表明 模拟 
与 实际 高 度 相符 '。 通 过 将 模拟 结果 与 2020 年 实 
际 土地 利用 数据 进行 对 比 ( 图 3)。 结 果 发 现 , Kappa 
系数 为 0.86, 模 拟 精 度 较 高 ,说明 PLUS 模 型 对 山西 
省 土地 利用 /土地 禾 盖 变化 模拟 是 有 效 的 ,可 用 于 模 
拟 山 西 省 未 来 土地 利用 变化 状况 。 


3 结果 与 分 析 


3.1 多 情景 下 土地 利用 需求 预测 

基准 情景 .城镇 化 情景 和 生态 保护 情景 3 种 情 
景 下 ( 表 7) ,林地 增长 面积 分 别 为 -4.67 km?、-4.74 
km ,23.94 km' ,其 中 ,林地 和 草地 的 面积 占 比 总 体 上 
保持 在 30% 左 右 ,耕地 面积 占 比 较 大 ,基本 维持 在 
35% 左 右 。 此 外 ,林地 内 部 转换 频率 较 高 , 相 比 于 灌 
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更 高 。 

基准 情景 是 对 研究 区 现 有 土地 利用 状况 的 顺 
延 ,在 这 种 情景 下 各 类 林地 的 面积 总 体 较 为 稳定 ， 
草地 为 有 林地 主要 转 人 对 象 。 耕 地 流失 面积 主要 
转向 城市 建设 用 地 ,表明 城镇 化 与 耕地 保护 的 矛盾 
已 经 成 为 山西 省 发 展 的 潜在 隐患 。 生 态 保 护 情景 
与 2020 年 实际 土地 利用 类 型 相 比 ,研究 区 内 林地 得 
到 很 好 的 保护 。 各 类 林地 面积 总 计 约 增加 23.94 
km2, 占 当前 总 林地 面积 的 5.38% ,有 林地 是 唯一 扩 
张 林 地 ,而 其 他 3 类 林地 都 有 小 幅 缩减 ,主要 原因 是 
林地 内 部 之 间 转 换 率 增 强 ,灌木 林 an PRO AS fi 
林地 向 有 林地 转 入 的 面积 分 别 为 1.74 km’, 10.28 
km^,0.37 km?。 城 镇 化 发 展 情景 是 基准 情景 的 强化 ， 
建设 用 地 增幅 较 大 ,与 2020 年 实际 相 比 ,增加 了 
10.79 km ,林地 除 灌木 林 外 总 体 上 出 现 了 不 同 程度 
的 退化 ( 表 4)。 
3.2 多 情景 下 土地 空间 分 布 模拟 

将 3 种 情景 下 山西 省 2030 年 土地 利用 面积 模 
拟 值 输入 到 PLUS 模 型 中 ,3 种 情景 下 林地 扩张 区 域 
主要 集中 在 研究 区 东南 部 长 治 盆地 附近 .中 西部 太 
原 盆 地 和 北部 大 同 盆地 。 有 林地 的 面积 扩张 范 于 
较 大 一 些 , 灌 木林 占 地 面积 总 体 上 保持 不 变 ,而 研 
究 区 西北 部 黑 茶 山 国家 自然 保护 区 和 东部 运城 分 
地 附近 玻 林地 扩张 变化 更 为 显著 ,其 他 林地 面积 小 
幅 降 低 。 总 体 来 看 ,东南 部 有 林地 呈现 大 面积 成 片 
扩张 ,中 西部 次 之 ,北部 林地 分 布 呈 碎片 化 特点 。 
受气 候 和 地 形 影 响 , 越 靠近 南部 地 区 ,林地 恢复 和 
扩张 的 可 能 性 越 大 (图 4)。 在 理想 状态 下 (生态 保 
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(a) 有 林地 (b) 灌木 林 
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(d) 其 他 林地 
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(c) IB PK 
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O 未 利用 地 


m 


图 2 2010 一 2020 年 各 土地 利用 类 型 适宜 性 概率 
Fig.2 Probability of suitability by land use type 2010-2020 


护 情景 ) ,通过 多 次 对 相关 参数 进行 调整 ,发 现 林 地 
总 面积 在 468.59 km 上 下 浮动 。 我 们 猜想 这 可 能 与 
山西 省 生态 环境 和 人 类 活动 等 因素 有 关 ,制约 了 林 
地 进一步 发 展 的 潜力 。 我 们 估 测 山西 省 2030 年 林 
地 总 体 增长 空间 约 为 23.94 km? , 占 当 前 林地 总 面积 
HY 5.38%. 
3.3 林地 扩张 典型 区 域 分 析 

生态 保护 情景 是 今后 山西 省 森林 修复 工程 的 
主要 参考 ,通过 提取 较为 明显 的 林地 扩张 区 域 (A、 
B.C AID) (AIS) ,并 将 其 与 其 他 情景 下 的 林地 变化 
进行 对 比 , 界 定 出 未 来 山西 省 林地 扩张 可 能 的 优先 
发 展区 (图 6) ,从 而 提高 山西 省 植树 造林 活动 效率 
和 固 碳 水 平 。 此 外 ,由 于 采取 较 高 分 辩 率 遥感 地 


图 , 故 对 于 林地 优先 发 展区 识别 可 以 精确 到 县 级 层 
面 ,这 也 是 本 文 主要 特点 。 

A 区 主要 位 于 人 研究 区 北部 ,包括 大 同市 和 朔 州 
市 ,该 区 域 林地 旦 碎片 化 分 布 ,从 林 分 变化 来 看 ,A 
区 域 总 体 表现 为 耕地 向 有 林地 的 转换 ,其 他 林地 转 
换 相对 稳定 。 该 地 区 有 林地 的 扩展 趋势 较为 明显 ， 
且 在 山 阴 县 出 现 了 大 片 的 有 林地 增长 。 其 他 县 域 
也 出 现 了 不 同 程度 的 林地 增长 ,但 相对 分 散 。 故 在 
未 来 10 a 中 该 区 应 提升 植树 造林 的 集约 化 ,效率 化 
水 平 。 

B 区 主要 分 布 在 研究 区 中 西部 ,涉及 忻州 市 . 太 
原市 和 吕梁 地 区 ,该 区 域 以 灌木 林 为 主 。 与 其 他 情 
景 相 比 , 岚 县 有 林地 出 现 了 成 片 增长 的 趋势 ,而 在 
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图 3 山西 省 2020 年 土地 利用 实际 结果 和 模拟 结 


Fig. 3 Land use simulation results for Shanxi Province 


交 城 县. 文 水 县 等 境内 更 适宜 灌木 林 生 长 。 在 生态 
保护 情景 下 , 交 城 县 边缘 地 带 出 现 了 小 面积 的 有 林 
地 退化 为 疏 林 地 ,我 们 猜想 可 能 是 在 该 情景 中 ,有 
林地 邻 域 权重 达到 峰值 ,出 现 林地 转换 的 可 能 性 更 
大 。 这 也 需要 引起 我 们 的 重视 ,未 来 该 地 区 也 许 会 
出 现 林 地 退化 的 可 能 性 。 

C 区 覆盖 了 晋城 市 南部 ,临汾 市 东北 部 以 及 长 
治 市 西北 部 ,结合 山西 省 地 形 图 可 知 ,该 区 域 主要 
位 于 太原 盆地 和 临汾 盆地 。 受 暖 湿 气 候 影 响 ,C 区 
是 山西 省 境内 林木 主要 集中 区 , 且 以 有 林地 为 主 。 
在 不 同情 景 下 ,C 区 各 类 林地 之 间 的 转换 均 呈 现 稳 
定 状 态 。 值 得 注意 的 是 位 于 祁 县 境内 的 琉 林地 变 
化 ,与 城镇 发 展 情景 相 比 ,生态 保护 情景 中 有 林地 
在 该 区 域 有 较 大 扩展 潜力 ,大 面积 的 牙 林 地 转变 为 


ChinaXiv 合 作 期 刊 


1 期 du ， 娜 等 :山西 省 林地 扩展 优先 区 识别 117 


有 林地 ,表明 良好 的 生态 保护 措施 能 够 推动 区 域内 
林木 扩展 。 上 述 区 域 在 未 来 林地 扩展 潜力 更 大 。 

D 区 大 部 分 位 于 晋城 市 东南 部 , 少 部 分 位 于 运 
城市 东南 部 。 该 地 区 也 是 阳 城 荞 河 猕猴 国家 级 自 
然 保 护 区 和 太 宽 河 国家 级 自然 保护 区 重点 区 域 ， 
此 总 体 生态 环境 基础 较 好 。 有 林地 主要 分 布 在 陵 
川 县 境内 , 且 不 在 同情 景 中 林地 之 间 的 内 部 转换 较 
为 稳定 ,我们 猜想 这 可 能 与 地 方 已 有 的 生态 环境 基 
础 有 关 。 而 平顺 县 境内 更 适合 草地 生长 ,在 不 同情 
景 中 该 地 区 草地 面积 一 直 较 为 稳定 ,而 林地 扩展 却 
更 为 缓慢 。 长 治 县 和 壶 关 县 境内 出 现 了 大 面积 的 
有 林地 增长 ,在 其 他 发 展 情景 中 ,这 两 个 县 主要 以 
下 林 地 为 主 , 在 生态 保护 情景 下 ,调整 了 林地 的 参 
数 指标 后 ,出现 了 大 面积 的 有 林地 增长 , 即 该 县 域 
下 林 地 和 有 林地 的 转换 更 为 明显 ,上 述 区 域 有 林地 
扩展 的 可 能 性 比较 大 。 


4 讨论 


41 关于 植树 造林 与 森林 修复 

森林 修复 仍然 是 当前 缓解 气候 变化 的 最 有 效 
战略 之 一 ,与 传统 的 “生态 修复 "不 同 , 基 于 生态 安 
全 格局 的 植树 造林 活动 更 强调 治理 和 修复 的 “主动 
适应 ”*。 和 干旱 地 区 受 自然 条 件 限制 ,植树 造林 活 
动 的 有 效 性 有 待 商检 ”。 因 此 ,如 何 提升 造林 的 
科学 性 和 “主动 适应 ”性 成 为 当下 山西 省 生态 保护 
的 重要 议题 。 从 结果 来 看 ,山西 省 森林 生态 系统 似 
乎 并 不 十 分 稳定 ,在 理想 的 状态 (生态 保护 情景 ) 
下 ,局 部 地 区 和 森林 边缘 已 经 出 现 了 和 森林 退化 的 趋 
势 ,如 果 不 调整 林地 领域 权重 值 , 这 种 情况 也 许 会 
更 严重 。 耕 地 仍然 是 山西 省 主要 的 土地 类 型 ,这 对 


表 7 多 情景 下 山西 省 土地 利用 面积 对 比 


Tab.7 Comparison of land use areas in Shanxi Province under multiple scenarios /km 
土地 利用 类 型 2020 年 实际 情景 2030 年 基准 情景 2030 年 生态 保护 情景 2030 年 城镇 化 情景 
有 林地 193.67 192.33 230 190.32 
灌木 林 167.94 164.97 166.2 168.16 
[SU 7421 74.21 63.93 73.09 
其 他 林地 8.83 8.47 8.46 8.34 
耕地 578.01 572.94 554.37 553.41 
草地 443.01 435.48 433.09 451.39 
水 域 15 12.02 12.43 12.18 
建设 用 地 86.01 93.34 85.55 96.87 
未 利用 地 1.05 0.95 0.68 0.96 
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图 4 2030 年 多 情景 下 林地 变化 模拟 


Fig.4 Simulation of forest land change in 2030 under multiple scenarios 
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图 5 生态 保护 情景 下 林地 扩张 
Fig. 5 Woodland expansion under the 


ecological protection scenario 


构造 稳定 的 森林 生态 系统 产生 胁迫 , 需 考 虑 耕地 红 
线 这 一 基本 底线 。 在 理想 状态 下 山西 省 未 来 植树 
造林 的 发 展 空 间 仅 有 5.38% , 而 现实 往往 比 理想 状 
态 更 糟糕 。 在 当前 的 气候 和 政策 背景 下 ,山西 省 植 
树 造 林 的 潜力 有 待 商 椎 ,至少 在 当下 不 适合 大 规模 
种 植 , 目 前 尚 能 解决 的 是 在 有 限 的 造林 空间 内 进 一 


步 实现 林地 优化 。 

从 结果 分 析 ,A 区 应 划 为 林地 扩展 一 般 区域 。 
林地 边界 形状 较为 复杂 ,破碎 化 程度 高 , 受 人 类 活 
动 的 频繁 干扰 ,可 能 存在 建设 用 地 扩张 的 胁迫 。 在 
新 一 轮 的 植树 造林 活动 中 需 综合 两 个 片区 的 人 口 、 
经 济 等 社会 因素 确定 各 区 域 林地 扩张 的 弹性 边界 ， 
并 加 强 对 汐 木 林 和 草地 的 种 植 。B 区 现 有 灌木 林 面 
积 较 大 ,虽然 存在 玻 林 地 向 有 林地 转换 的 潜力 ,但 
是 从 因地制宜 的 角度 分 析 ,灌木 林 可 能 更 适宜 该 地 
区 的 气候 条 件 。C 区 和 D 区 应 划分 为 林地 扩展 极 重 
BED o ML GR UR SEC BONE , 碳 储 量 发 展 空间 较 大 。 
在 温带 和 亚热带 森林 中 ,树种 的 多 样 性 对 碳 储量 的 
增加 产生 积极 影响 ”“。 因 此 对 C 区 和 D 区 的 林地 扩 
展 应 注重 树种 多 样 性 的 提升 。 

4.2 研究 不 足 

最 后 需要 指出 ,本 研究 今后 至 少 需要 在 以 下 两 
个 方面 继续 推进 : (1) 模型 方面 ,PLUS 模型 虽然 可 
以 模拟 多 种 情景 下 土地 利用 变化 情况 ,但 是 基于 马 
尔 科 夫 链 模拟 多 种 发 展 情景 一 定 程度 上 具有 主观 
性 ,对 人 研究 区 经 济 发 展 的 实际 情况 欠 考 虑 。 虽 然 土 
地 利用 变化 特征 也 有 较为 明显 的 马尔 科 夫 特征 ,但 
模型 事实 上 并 未 考虑 社会 经 济 发 展 以 及 不 同 生产 
部 门 对 土地 的 需求 。 而 且 若是 不 同时 段 内 发 展 十 
地 变化 速度 有 较 大 差异 ,预测 的 结果 也 有 很 大 的 差 
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生态 保护 情景 
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图 6 林地 优先 发 展区 识别 


Fig.6 Identification of priority woodland development areas 


异 。 在 诸多 研究 中 ,有 关 学 者 采用 系统 动力 学 模型 
来 推测 不 同 发 展 情景 下 土地 利用 场景 ” ,更 具备 严 
间 性 ,这 也 为 本 研究 的 进一步 延 申 提供 了 参考 ,(2) 
数据 方面 ,关于 识别 框架 的 指标 选取 ,影响 森林 扩 
张 的 因素 有 很 多 ,例如 气候 条 件 .土壤 肥力 等 “”， 
但 目前 尚未 有 可 用 的 高 精度 空间 数据 进行 分 析 。 


5 结论 


本 研究 讨论 了 在 所 有 造林 的 规划 阶段 遇 到 的 


两 个 关键 性 问题 :一 是 何 处 是 造林 适宜 区 ”二 是 
如 何 集 中 ?” 本 研究 构建 了 针对 山西 省 林地 优先 发 
展区 识别 框架 对 上 述 问 题 进行 分 析 , 借 助 PLUS 模 
型 和 马尔 科 夫 链 系统 模拟 了 2030 年 山西 省 不 同情 
景 下 的 土地 利用 变化 情况 ,并 界定 未 来 林地 潜在 优 
先 区 。 研 究 得 出 以 下 结论 : COD 在 当前 气候 和 政策 
条 件 下 ,山西 省 2030 年 林地 增长 空间 非常 有 限 , 仅 
为 5.38%。 主 要 宜 林 地 区 域 集中 在 东南 部 ,以 有 林 
地 扩张 为 主 ;其 次 为 中 部 地 区 ,以 灌木 林 扩 张 为 主 ， 
且 林 地 边缘 面临 有 林地 退化 的 风险 ;最 后 为 北部 ， 
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基于 林地 破碎 化 严重 的 现状 ,未 来 修复 的 导向 更 偏 


向 于 草地 扩张 。(2) 林地 内 部 转换 频率 较 高 ,各 类 林 
地 按照 扩展 潜力 高 低 划 分 ,依次 为 :有 林地 > 灌木 
林 > 政 林 地 > 其 他 林地 。(3) 从 当前 山西 省 的 发 展 趋 
势 来 看 ,土地 利用 的 方式 是 非 可 持续 的 ,我 们 对 山 
西 省 大 规模 植树 造林 的 方式 保持 谨慎 态度 。 且 在 
耕地 保护 红线 的 限制 下 ,认为 生态 保护 情景 的 发 展 
模式 在 山西 省 实施 起 来 可 能 会 存在 困难 。 
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Identification of priority areas for forest land expansion in Shanxi Province 


MENG Na, ZHANG Ying 


(School of Economics and Management, Beijing Forestry University, Beijing 100083, China) 


Abstract: Over the past 70 years, large-scale afforestation activities have been conducted in arid and semi-arid 
areas in China with remarkable success, but localized afforestation has failed due to regional human activities and 
climatic influences. Although the potential for afforestation in arid areas is recognized, it is unknown in which 
micro-region afforestation should be concentrated. Taking Shanxi Province, a typical region of the Loess Plateau, 
as an example, this study constructs a framework for the identification of forestable land and uses the PLUS and 
Markov models to simulate land use changes under different scenarios in 2030, analyzing woodland expansion, 
trends in internal changes in forest stands, and spatial distribution. The results show that: The current space for 
woodland growth in Shanxi Province is 5.38%, and the southeast is the main potential growth area; there is a 
possibility of woodland degradation at the edges of woodlands in the central and western parts of the province, 
while the degree of woodland fragmentation in the north is higher; the frequency of intra- forest conversion is 
higher; in order of expansion potential, there are woodland > shrubland > open woodland > other woodlands. In 
this context, a cautious attitude is taken toward large-scale afforestation in Shanxi Province. This study can 
provide a reference for effective afforestation management and enhancement of forest carbon sequestration levels 
in the Shanxi Province. 

Keywords: PLUS model; Markov model; afforestation; scenario simulation; priority zone identification; arid 


and semi-arid areas 


